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We theoretically discuss the mechanism for the peculiar nonlinear conduction in quasi-two-dimensional or- 
ganic conductors 6I-(BEDT-TTF)2X [BEDT-TTF=bis(ethylenedithio)tetrathiafulvalene] through the melt- 
ing of stripe-type charge order. An extended Peierls-Hubbard model attached to metallic electrodes is in- 
vestigated by a nonequilibrium Green's function technique. A novel current-voltage characteristic appears 
in a coexistent state of stripe-type and nonstripe 3-fold charge orders, where the applied bias melts mainly 
the stripe-type charge order through the reduction of lattice distortion, whereas the 3- fold charge order sur- 
vives. These contrastive responses of the two different charge orders are consistent with the experimental 
observations. 
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Nonlinear conduction in low-dimensional electron sys- 
tems has been of great interest from the viewpoint of 
fundamental nonequilibrium physics and possible appli- 
cations to electronic devices. A well-known example is a 
sliding of density waves in quasi-one-dimensional materi- 
als, where a nesting of the Fermi surface is responsible for 
their ground states.^) In strongly correlated systems such 
as Mott insulators^-' and charge-ordered states of tran- 
sition metal oxides,'^' dielectric breakdown phenomena 
have been observed. In one dimension, a breakdown of 
Mott insulators by the Landau-Zener tunneling mecha- 
nism has been proposed theoretically^-* and the relevance 
to experimental findings has been discussed so far. 

The observations of giant nonlinear conduction and 
spontaneous current oscillation in the organic com- 
pounds 6l-(BEDT-TTF)2CsM(SCN)4'5'S) (M=Co and 
Zn) have renewed interest since their mechanism and the 
electric-field-induced behaviors seem to differ in many re- 
spects from those in the above materials. The family of 
organic conductors 6'-(ET)2X (ET is the abbreviation of 
BEDT-TTF) is known to exhibit charge order (CO).'^'®) 
It has a quasi-two-dimensional structure, where ET 
' molecules form a triangular lattice in each conduction 
layer [Fig. l(a)]^^ with electron density at 3/4 filling (one 
hole per two ET molecules). In 6'-(ET)2RbZn(SCN)4, a 
metal-insulator transition with a structural distortion oc- 
curs at Tco = 190 K.^^ A stripe-type arrangement of lo- 
calized charges indexed by wave number — (0, 0, 1/2), 
which is called a horizontal CO [Fig. 1(b)], emerges below 
Tco-^°"^^^ In 6l-(ET)2CsM(SCN)4, on the other hand, 
coexistence of two different COs, the horizontal CO and 
a nonstripe CO indexed by qi = {2/3, k, 1/3), has been 
observed by X-ray experiments,^'^' ^^"^ although this com- 
pound does not show any long-range order. Several ex- 
periments suggest that the peculiar nonlinear conduc- 
tion results from suppression of the horizontal CO by an 
electric field without destroying the nonstripe CO.^'^^^ 
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A similar nonlinear conduction and coexistence of two 
kinds of COs have also been observed in a rapidly cooled 
RbZn salt.^^' Although nonlinear conduction has been 
found in other compounds,^®' ^^-^ this unique feature of 
multiple charge modulation seems essential for the non- 
linearity in ^-type compounds. 

Theoretically, the CO phenomenon in 6'-(ET)2X has 
been investigated from various aspects. ^'^'^^^ The charge 
disproportionation results mainly from the long-range 
nature of the Coulomb interaction, whereas electron- 
phonon (e-ph) couplings are also important. In partic- 
ular, the lattice distortion in the RbZn salt considerably 
stabilizes the horizontal CO.^^"^"*' Compared with our 
knowledge on the ground states, nonequilibrium states 
induced by an external field have been poorly under- 
stood. Recently, Mori and coworkers have shown that a 
phenomenological equation can reproduce the observed 
nonlinear current-voltage characteristics in some com- 
pounds.^^' However, the origin of nonlinear conduc- 
tion is still unclear, so a microscopic theory is highly 
desirable. 

With these in mind, we investigate nonequilibrium 
steady states of 0-(ET)2X under applied bias voltages, 
using a model that takes account of both the long- 
range Coulomb interactions and e-ph couplings. ^^"^''^ 
The model describes competition among various COs, 
including the horizontal CO and a so-called 3-fold CO 
[Fig. 1(c)], the latter has a nonstripe charge pattern^^^ 
and can be related to the CO with qi in the CsM salt. 
The state with the horizontal CO is insulating, whereas 
that with the 3-fold CO is metallic. ^^^^ We show that 
when these two states coexist, the bias voltage melts the 
horizontal CO and largely alters the conduction behav- 
ior. The lattice distortion, which induces the horizontal 
charge modulation, has a key role in determining whether 
the system becomes resistive or conductive. 

We consider the extended Peierls-Hubbard model^^"^^^ 
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Fig. 1. (Color online), (a) Schematic view of the high- 
temperature structure for 6-(ET)2X. Ellipses represent ET 
molecules. The gray rectangle shows a unit cell in which sites are 
labeled A and B. (b) Horizontal CO and (c) 3-fold CO, where the 
solid (open or shaded) ellipses denote hole-rich (-poor) molecules. 
In (b), the lattice distortion along the horizontal stripe is shown 
by the thick and broken lines. 

written as 

+ C/^(n,t-iVe/2)(n,^-A^e/2) 

i 

where (ij) represents the summation over pairs of neigh- 
boring sites, cl^{cicr) denotes the creation (annihilation) 
operator for an electron with spin a at the ith site, 
n-ia = cl^Cia, Tii = Tii^ + Tii^, and the averaged electron 
density Ng = 1.5. tij denotes the transfer integrals and 
U the on-site repulsion. For the intersite Coulomb inter- 
actions Vij, we consider up to third-neighbor pairs of 
sites, the summation over which is represented by {{ij))- 
The e-ph coupling constant, the lattice displacement and 
the elastic constant are denoted by aij, Uij, and Kij, 
respectively. We introduce new variables as yi j = ai jUi j 
and s„ =a?,^./i^.,,. 22-24) 

The structure of ^-(ET)2X in the high-temperature 
metallic phase is shown in Fig. 1(a). There are two trans- 
fer integrals tc and tp on the vertical and diagonal bonds, 
respectively. We set tc = —0.04 (eV) and tp = 0.1 in 
the following. On the vertical (diagonal) bonds, we de- 
fine the nearest-neighbor interaction Vc (Vp). For the 
second- and third-neighbor interactions, we write them 
as Vi,j — Vir/vij. Here, rij is the distance between the 
ith. and jth sites. The horizontal CO and the 3-fold CO 
are schematically shown in Figs. 1(b) and 1(c), respec- 
tively. For e-ph couplings, we consider a lattice distortion 
caused by the molecular rotation, which is crucial for sta- 
bilizing the horizontal 00.^2-24) xhis is because it gives 
homogeneous modulation in the transfer integrals along 
with a horizontal stripe. We denote the modulation as j/^ 
and write the corresponding e-ph coupling as s^. These 
variables are defined only on the bonds that are depicted 
by the thick and broken lines in Fig. 1(b). We assume 
that is independent of bond index so that the transfer 
integrals on the hole-rich (-poor) stripe are written as 
tp + y<i> {tp — y^) with > 0. The elastic energy for each 
distorted bond is given by y'^/{2s^). For simplicity, we 
do not take account of other types of modulation. 22-24) 
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Fig. 2. (Color online). Schematic picture of the model. The left 
and right electrodes are attached to the central part where the 
horizontal and 3-fold COs coexist when the bias voltage is absent. 
A horizontal charge modulation induced by the lattice distortion 
on the 3-fold charge pattern is represented by the gray ellipses on 
the thick bonds. The x-axis (y-axis) is along (perpendicular to) the 
conduction direction. 

We describe the steady states under applied bias 
voltages by the nonequilibrium Green's function 
method. 2'^' 28) gj^gwn in Fig. 2, we attach the left and 
right (a = L,R) semi-infinite metallic electrodes to the 
central part that is described by eq. (1). A coexistent CO, 
which will be discussed later, is realized in the absence of 
the bias voltage. For the ith site [i = {ix, iy)] in the cen- 
tral part, we define its coordinates as {ix, iy) if ix is odd, 
and {ix, iy — 1/2) if ix is even. The numbers of sites are 
denoted by Lx and Ly (1 < ij; < Lx and 1 < iy < Ly). 
We assume that electrons in the leads are noninteract- 
ing and that they move only in the x-direction; for sim- 
plicity, the electrodes are one-dimensional. The effects 
of the leads on the central part are incorporated into 
the self-energies. In the wide-band limit, the retarded 
self-energies are independent of energy and written as 
(E^)y = -^jaJ2i^ ^lic^jic^ where % is the Kronecker 
delta, iL {in) denotes the site that is connected with 
the left (right) electrode, and the coupling constant 
between the central part and the electrode a. 

We use the Hartree-Fock approximation for the U, 
Vp, and Vc terms in eq. (1). For the second- and third- 
neighbor interactions, we employ the Hartree approxima- 
tion. These interactions contribute to a redistribution of 
charges near the electrodes, 2^) whereas the charge dispro- 
portionation is mainly caused by Vp and Vc- We consider 
the case of Vp, Vc ^ Vr- The periodic boundary condi- 
tion is adopted along the y-axis. In the mean-field Hamil- 
tonian, we introduce a scalar potential tp that is defined 
by the Hartree terms for the intersite Coulomb interac- 
tions as ■>p{ix, iy) = ^i,3{{'^j) - ^e) + aix + b, which 
is equivalent to the Poisson equation. 2''') The slope of the 
potential, aix, describes the electric field in the central 
part. 2^) The constants a and b are so determined that -0 
satisfies the boundary conditions, '4'{1, iy) = V/2 
and ■^J2iy''P{^x,iy) = —V/2, when the bias voltage 
V is applied to the system. We assume that the work- 
function differences at the interfaces are absent. 

The steady states under the applied bias are obtained 
as a self-consistent solution for the mean fields that 
are calculated by the method described in ref. 27. The 
electron density {riia-) is calculated by decomposing it 
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Fig. 3. (Color online), (a) Energies of the horizontal (hr) and 
coexistent (coex) COs with Lx = Ly = 18, U = 0.6, Vc/U = 0.35, 
= 0.1, Vir = 0.02, and 7^ = 7l = 0, relative to that of the 
3-fold CO, as a function of Vp/Vc- (b) The lattice distortions y^, 
for the horizontal and coexistent COs. Density of states for (c) 
the 3-fold, (d) horizontal, and (e) coexistent COs in the case of 
Vp/Vc = 1. 



Fig. 4. (Color online), (a) Electric current J, (b) lattice distor- 
tion and (c) charge structure factors 5c (qhr) s-i^d ■^cCqaf) as 
a function of V, where q^r and q3f are the wave vectors for the 
horizontal (hr) and 3-fold COs, respectively. In (a) and (b), results 
with Lx = Ly = 24 are also shown. In (c), the qy dependence of 
Sc{0,qy) is shown in the inset, (d) Density of states for V = 0.08 
and 0.21. 



into the "equilibrium" and "nonequilibrium" parts as 
{nicr} — n'^i'^ + X]q '^'^icT-^'^'' The same decomposition is 
used to obtain {c\„Cjcr). We adjust the chemical poten- 
tial /i(7 — (/iL-l-/ifl)/2^^) such that the electron density of 
the central part is fixed at 3/4 filling. Here, /i/, and 
are the left and right chemical potentials, respectively. 
For finite V, we set = fi(j + V/2 and = fJ-c — V/2. 
The lattice distortion is determined as in the equi- 
librium case.^^-* The current J is obtained by using Snf^ 

as J = 7fl ^rif^a - IL Y.^^a ^^f^J'^^ ^hcrC WC SCt 

e = ?i = 1. In the following, we use U — 0.6, VcjU = 0.35, 
S0 = 0.1, and Vir — 0.02. The size of the central part is 
Lx — Ly ^ 18 unless otherwise noted. 

First, we consider the equilibrium case where the CO 
system is isolated from the electrodes (7^ = 7l = 0). 
The ground-state energies per site of three mean-field 
solutions as a function of Vp/Vc are shown in Fig. 3(a), 
where the energy of the 3-fold CO is set at zero. The hor- 
izontal and coexistent COs have a finite lattice distortion 
t/0, as shown in Fig. 3(b). For Vp/Vc < 0.7, the horizon- 
tal CO is the most stable, whereas the 3-fold CO has a 
lower energy than the horizontal CO for Vp/Vc > 0.7. 
This results from charge frustration on the triangular 
lattice.25'26) The 3-fold CO is further stabilized by co- 
existing with the ?/0-induced horizontal CO (Fig. 2) and 
becomes the ground state near Vp/Vc = 1-^^^ For each 
CO pattern, we show the density of states D{E) in Figs. 
3(c)-3(e) for = 1, where we used a broadening fac- 

tor of ry = 0.01. A finite D{E) at the Fermi level exists 
for the 3-fold CO since this state is metallic, whereas 
the horizontal CO has an energy gap. In the coexistent 
CO, D{E) at E = is suppressed compared with that 
for the 3-fold CO, which is due to the horizontal CO. 
However, we note that the state has no energy gap in 
the thermodynamic limit. 



Next, we discuss the results obtained with finite bias 
V. We set Vp/Vc = 1 and 7^=7^ = 0.03, where the 
ground state without the bias is the coexistent CO. The 
current-voltage characteristics and the distortion are 
shown in Figs. 4(a) and 4(b), respectively. In these fig- 
ures, we show the results with L^ = Ly = 24 for com- 
parison. The behaviors of J and are qualitatively the 
same as those with L^ = Ly = 18. With increasing V, the 
current gradually increases and abruptly becomes large 
at Vcr ^ 0.17. For small V and L^. — Ly — 18, J has 
stepwise structures owing to the finite-size effect. As 
shown in Fig. 4(b), j/^ is almost unchanged for V < Vcr, 
although it shows a gradual decrease with increasing V. 
At V = Vcr, 2/0 steeply decreases, which is directly re- 
lated to the reduction of the horizontal CO. 

The effects of bias voltages on the horizontal and 3- fold 
components of the charge distribution are obtained by 
calculating the charge structure factor, which is defined 
as 

^c(q) ^^Y.^{n^^n,A) + {n^.Bn.B))e'^^^--'^^l (2) 

Here, we use the unit cell shown in Fig. 1(a), which is la- 
beled by fx, V. A and B are indices for sites inside the unit 
cell. The position vector for the /i-th (i^-th) unit cell is 
denoted by R^, (R-y) and Ng ~ L^xLy. The wave vectors 
that correspond to the horizontal and 3-fold components 
are qhr ~ (0, n) and qaf = (0, 27r/3), respectively. In Fig. 
4(c), we show S'c(qhr) and S'c(q3f) as a function of V, 
and the qy dependence of Sc{0,qy) for y = and 0.21 
in the inset. In the V = case, Sc{0, Qy) has two peaks 
at Qhr and qsf since the two COs coexist. Because the 3- 
fold charge modulation is larger than the horizontal one, 
we have S'c(q3f) > 5*0 (qhr)- For V < Vcr, the values of 
5'c(qhr) and S'c(q3f) are almost unchanged, which indi- 
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cates that the coexistent CO is robust against the apphed 
bias. For V > Vcr, both S'c(qhr) and 5'c(q3f ) decrease and 
the charge distribution is largely modified. In particu- 
lar, the horizontal component is drastically weakened. 
For V = 0.21, the qhr peak in Sc{0,qy) disappears. Cor- 
respondingly, becomes very small for large V. How- 
ever, the 3-fold component survives even in the region 

V > Vcr- It is noteworthy that J increases without de- 
stroying the 3-fold CO. We show the density of states for 

V = 0.08 and 0.21 in Fig. 4(d). For V = 0.08, D{E) is 
qualitatively the same as that in the V = case shown in 
Fig. 3(e). Since the coexistent CO has no energy gap,^^^ 
a small current can flow even for V < Vcr- For V = 0.21, 
a large D{E) appears at E = jic since the horizontal 
charge modulation is suppressed. The change in the con- 
duction behavior is triggered by the reduction in y^. If 

is large, the horizontal CO persists so the system is resis- 
tive, whereas if y^ decreases, only the 3-fold CO remains 
so the system becomes conductive. 

Let us disciiss the relevance to experimental results on 
^-(ET)2X. Om results basically reproduce the X-ray re- 
sults^' that indicate that the nonlinear conduction 
is caused by the melting of the horizontal CO whereas 
the nonstripe CO remains. In the CsZn salt, the resis- 
tivity begins to increase at around 50 K.^^ This cor- 
responds to the growth of X-ray intensity for the hor- 
izontal CO, whereas that for the nonstripe CO shows 
only slight temperature dependence. "'^^^ Theoretically, 
such contrastive temperature dependences of stripe-type 
and 3-fold-type charge fluctuations have been shown by 
the random phase approximation,^^) where the former 
comes from the Fermi-surface instability assisted by e- 
ph couplings, whereas the latter is due to the wave- 
vector dependence of the Fourier transform of the inter- 
site Coulomb interaction. These facts also suggest that 
the horizontal CO induced by the lattice distortion is 
directly related to the resistive behavior, which is consis- 
tent with our results. However, there are some issues that 
remain to be clarifled. Experimentally, the nonlinear con- 
duction appears in a state with no long-range CO, which 
is in contrast to the results of our mean- field calculations. 
Therefore, at present, it is difficult to compare the results 
quantitatively. In fact, the electric fields required for non- 
linearity (~1 V/cm for the CsZn salt and ~10 V/cm for 
the rapidly cooled RbZn salt) are much smaller than that 
in the present study. For a quantitative comparison, the 
effects of quantum fluctuations must be taken into ac- 
count. The finite-size effects^^^ as well as the thermal 
fluctuations will affect the values of the threshold volt- 
age. Recently, an inhomogeneous state of competing COs 
has been suggested and a possible relation to the nonlin- 
ear conduction has been discussed.'"'' ^'^"^^^ Although the 
present study is based on a uniform CO, we speculate 
that in a spatially nonuniform state, only domains of the 
horizontal CO arc suppressed by an electric field, which 
resiilts in nonlinear current-voltage characteristics sim- 
ilar to our results. The origin of the current oscillation 
in (?-(ET)2X'''' occurring with the nonlinear conduction 
is still unclear. However, Suko et al.^^ have recently sug- 
gested that the oscillation is due to a current-induced 



modulation of the lattice distortion, which may be re- 
lated to our results. 

In summary, we have investigated the mechanism of 
nonlinear conduction in 6'-(ET)2X. In the coexistent 
state of horizontal and 3-fold COs, the bias voltage weak- 
ens the lattice distortion and melts the horizontal CO. 
The metallic 3-fold CO remains even after the disap- 
pearance of the horizontal CO, which leads to selective 
melting of the latter CO. We have shown that these dif- 
ferent responses depending on the spatial patterns of the 
two COs bring about the novel nonlinear conduction. 
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